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Abstract: 



O ■ Future hadron coUider experiments are expected to record large to huge samples 
of ti events. The analysis of these data with respect to ti spin-spin correlations 

^ , requires precise predictions of the production of these quark pairs in a general spin 

p • configuration. Therefore we have computed, at next-to-leading order (NLO) in the 

p • QCD coupling, the spin density matrices describing ti production by quark antiquark 

D . annihilation, which is the dominant production process at the Tevatron. Moreover 

^ ' we have computed the strength of the ti spin correlation at NLO, using various spin 



quantization axes. 
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Pairs of top and antitop quarks will be produced copiously with the upgraded Fer- 
milab Tevatron collider and in even larger numbers with the Large Hadron Collider 
(LHC) at CERN, allowing for detailed investigations of the properties of these quarks 
in the future. For this aim the study of top spin phenomena will play an important 
role. This is because the top quark, as compared with lighter quarks, is unique in 
that it is sufficiently short-lived to prevent hadronization effects from diluting the 
spin-polarisations and spin-spin correlations that were imprinted upon the t and t 
quarks by their production mechanism. Hence spin-polarization and spin-correlation 
phenomena will provide valuable information about the interactions of top quarks. 
An attempt to detect these spin correlations in a small tt dilepton sample collected 
at the Tevatron was recently reported by the DO collaboration |T[]. 

Needless to say, in order to interpret future data, the predictions of these spin 
phenomena should be as precise as possible. Within the Standard Model (SM) top 
antitop pair production at hadron colliders is determined by QCD. So far the next- 
to-leading order (NLO) QCD corrections are known for the spin-averaged differential 
ti cross section 0, |^, ^, |^, while spin correlations were analysed only to leading- 
order P, 1^, H, ^ [l^! mi in In this letter we report on progress in providing a 



complete description of these spin effects within QCD at NLO: We have computed, 
as a first step in this direction, the ti spin density matrices for the parton reactions 
qq — > tt,tig to order af. We have also calculated, for these reactions, the degree of 
the ti spin correlation at NLO for different t and i spin quantization axes. 

At the Tevatron qq annihilation is the dominant process for producing top quark 
pairs while at the LHC ti production is mainly due to gluon gluon fusion. The ti 
spin correlations and spin polarisations can be inferred from appropriate angular 
correlations and distributions of the t and i decay products. In the SM the main 
top decay modes are t — > bW — * bqq', biue. Among these final states the charged 
leptons, or the jets from quarks of weak isospin -1/2 originating from W decay, are 
most sensitive to the polarisation of the top quarks. Hence to leading order in the 
QCD coupling one has to treat the two parton reactions 

gg,qq-^ti^bb + 4f, (1) 

where f = q,i, z/£. The amplitudes for these two- to six-body processes, with inter- 



mediate top quarks of non-zero total width, were given first in |12|. The computation 
of the NLO QCD corrections to these amplitudes is quite a demanding task In 
view of the fact that the total width of the top quark is much smaller than its 
mass, Tt/rrit = 0{1%), one may however analyse these reactions using the so-called 
narrow width or leading pole approximation |TB|, In our case this approach con- 
sists, for a given process, in an expansion of the amplitude around the poles of the 
unstable t and i quarks, which amounts to an expansion in powers of Tt/rrit. Only 
the leading term of this expansion, i.e., the residue of the double poles is considered 



^At NLO, apart from the gluon radiation reactions, the processes g + q{q) tt + q{q) 
bb + Af + q{q) must also be taken into account. 
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here. Then the radiative corrections to (|l|) can be classified into so-called factoris- 
able and non-factorisable corrections. The non-factorisable NLO QCD corrections 



were calculated in |15]. For the factorisable corrections the square of the complete 



matrix element Ai^'^^ is of the form 

I ^(^f (X Tr [pR^'^p] = Pa'aRi'l,,f,p>P^'(^. (2) 

Here i?*-^^ denotes the density matrix for the production of on-shell ti pairs, the 
label A indicates the process, and p, p are the density matrices describing the decay 
of polarised t and i quarks, respectively, into specific final states. The subscripts in 

denote the t, t spin indices. Note that both the production and decay density 
matrices are gauge invariant. 

The one-loop QCD corrections to the semileptonic decays of polarised top quarks 
and to t ^ W + b can be extracted from the results of |TB| and |TB|, respec- 



tively. In the following we describe our computation of the density matrices for tt 
production by qq annihilation. At NLO we have to consider the reactions 

q{pi)+q{p2)^t{h) + t{k2), (3) 

and 

qipi) + q{p2) ^ t{h) + t{k2) + gih). (4) 

We define the production density matrix for the process (|^) in terms of its transition 
matrix element as follows: (for ease of notation we omit the label A in the following) 

Raa',pp' = ^ ^ {taip\T\qq) {qq\T^Udp,) , (5) 

colors 
initial spins 

where the factor Nqq = {2Ncf = 36 averages over the spins and colors of the initial 
qq pair. The matrix structure of R is 



where cr' are the Pauli matrices. Using rotational invariance the 'structure functions' 
Bi, Bi and Cij can be further decomposed. The function A, which determines the 
it cross section, is known to next-to-leading order in from the work of [Q, ^. 
Because of parity (P) invariance the vectors B, B can have, within QCD, only a 
component normal to the scattering plane. This component, which amounts to a 
normal polarisation of the t and i quarks, is induced by the absorptive part of the 
scattering amplitude, and it was computed for qq and gg initial states in |]19|, ^ 
to order a^. The normal polarisation is quite small, both for ti production at the 
Tevatron and at the LHC. Parity and CP invariance of QCD dictates that the 
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functions Cij, which encode the correlation between the t and t spins, have the 
structure [Blf 



Cij = CiSij + C2PiiPij + cshikij + c^ikiiPij + puhj), (7) 

where pi and ki are the directions of flight of the initial quark and of the t quark, 
respectively, in the parton cm. frame. The production density matrix for the 
reaction (||) can be defined and decomposed in an analogous fashion. 

To Born approximation the functions were given, e.g., in 0. In order to 
determine these functions to order we first computed the one-loop diagrams 
that contribute to (^. Dimensional regularization was employed to treat both the 
ultraviolet and the infrared and collinear singularities which appear in the diagrams. 
The ultraviolet singularities were removed by using the MS prescription for the QCD 
coupling as and the on-shell definition of the top mass m^. The initial quarks are 
taken to be massless. After renormalisation the density matrix for the ti final state 
still contains single and double poles in e = (4 — D)/2 due to soft and collinear 
singularities. These poles are cancelled after including the contributions of the 
reaction @ and mass factorization. For the latter we used the MS factorization 
scheme. We avoided the computation of the exact density matrix for the reaction 

in D dimensions by employing a simple version of the phase-space slicing method 
[ P^ , We divided the phase space into four regions, namely the region where 
the gluon is soft, the two regions where the gluon is collinear to one of the initial 
state massless quarks (but not soft), and the complement of these three regions, 
where all partons are 'resolved'. This decomposition can be performed using a 
single dimensionless cut parameter Xmin- For example, the soft region is defined 
by the requirement that the scaled gluon energy in the cm. system Xg = 2Eg/ y/s 
is smaller than Xmin. In the soft region we used the eikonal approximation of the 
matrix element for reaction (^) and the soft limit of the phase space measure. The 
integration over the gluon momentum can then be carried out analytically in D 
dimensions. The two collinear regions are defined by (cos^^g > (1 — Xmin) and 
Xg > Xmin) and {cos 9qg < (—1 -|- Xrain) and Xg > Xmin), respectively, where 9qg is 
the angle between the gluon and the quark in the qq cm. frame. In these regions 
we used the collinear approximations for both the squared matrix element and the 
phase space in D dimensions. Finally, the exact spin density matrix for reaction (^) 
in four space-time dimensions was used in the resolved region, where all necessary 
phase space integrations can be carried out numerically. By construction, all four 
individual contributions depend logarithmically on the slicing parameter Xmin, but 
in the sum only a residual linear dependence on Xmin remains, which is due to the 
approximations made in the soft and collinear regions. By varying Xmm between 
10~^ and 10~^ we checked that for Xmin < 10"'^ this residual dependence is smaller 
than our numerical error (which is less than a permill for all results discussed below). 

After mass factorisation we are left with finite density matrices for the tt and the 
ti + hard gluon final states. As a check of our calculation we first compute the total 
cross section for qq tt + X at NLO. If one identifies the MS renormalisation scale 
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/i with the mass factorisation scale hf and neglects all quark masses except for rrit, 
then one can express the cross section in terms of dimensionless scaling functions 



'"'t 



where s is the parton cm. energy squared and rj = s/Amf — 1. We have compared 
our result for function of t] with those of p|, |^ and found perfect agreement. 

We now consider the following set of spin-correlation observables: 

(9) 

(10) 

(11) 
(12) 

(13) 

where Si,S2 are the t and t spin operators, respectively. The factor of 4 is conven- 
tional. With this normalization, the expectation value of Oi is equal to 1 at the 
Born level. The expectation values of the observables (p^Of), (p!l|), ([T^), and (|T3|) 
determine the correlation of different t, i spin projections. Eq. ([T0|) corresponds to 
a correlation of the t and i spins in the helicity basis, while (0) correlates the spins 
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projected along the beam line in the parton c.m.s. The 'beam-line basis' used in (12) 
was defined in and refers to spin axes being parallel to the antiquark direction in 
the t rest frame P2 and to the quark direction in the t rest frame p^*, respectively. 
The spin axes di 2 in (|l^) correspond to the so-called 'optimal basis' p4|, |10| to be 
discussed below. 

For quark-antiquark annihilation it turns out that the spin correlation (|TTD [0, |ri[| 
and the correlation in the beam-line basis (O) [H are stronger than the correlation in 



the helicity basis. A spin-quantization axis was constructed in ||2J, |I0[ with respect 
to which the t and i spins are 100% correlated to leading order in the QCD coupling, 
for all energies and scattering angles. In terms of the structure functions of @ this 
means that the 'optimal' spin axis d fulfills the condition 

diCjd, = A. (14) 

The existence of a solution of Eq. (|T^ is a special property of the leading order spin 



density matrix for the reaction qq ^ tt. One finds |23, ITO 



pi + (1 -7i)(pi -k^ki ^^^^ 



l-(pi-ki)2(l-7?) 
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where 71 = Ei/rrit. The construction of this axis exphcitly uses the leading order 
result for the spin density matrix, and different generalizations to higher orders are 
possible. We use in (O) as spin axes: 



di = d, 

-pi + (1 - 72)(pi ■ k2)k2 



(16) 



l-(pi -1.2)2(1-71 



2^ 



where 72 = E2/mt. For the 2 to 2 process qq ti, d2 = di = d. (A different 
generalization of d to higher orders was used in p5| .) 

The expectation value of a spin-correlation observable O at parton level can be 
written at next-to-leading order in analogy to (||) as follows: 

(O),, = gll'^iv) + ^^c^sigSiv) + -gSiv) Hf^l/m',)). (17) 

Note that these quantities depend exphcitly on the factorization scale ftp, but only 
implicitly (through a^) on the renormalization scale /x. This is because a factor 
drops out in the expectation values, and hence the Born result is of order a^. 

Our results for the functions gfq{rj), g^qq{r]) and g^qq{rj) are shown for the five 
observables (PD-(p!3|) in Figs. 1-5. In each figure, the dotted line is the Born result 
gfq{ri), the full line shows the function g^qqirj), and the dashed line is g^qq{rj). A 
general feature of all results is that the QCD corrections are very small for values 
of ?7 < 1. For larger values of r], the functions g^qq{ri) depart significantly from 

zero. Also, the functions gq]i{'ri) become nonzero, with less dramatic growth as 
1] ^ (yo and with an opposite sign as compared to g'fq{ri)- The phenomenological 
implications of these features for spin correlations at the Tevatron and the LHC will 
be studied in detail in a future work. Here we merely note that the substantial QCD 
corrections for large 77 will be damped by the parton distribution functions, which 
decrease rapidly with 77. Moreover, at Tevatron energies values of r] above ~ 30 are 
kinematically excluded. 

To summarize: We have computed the spin density matrices describing ti produc- 
tion by qq annihilation to order a^. Further we have evaluated the scaling functions 
encoding the QCD corrections to spin correlations, using a number of different spin 
quantization axes. This work provides a building block, which was missing so far, 
towards a complete description of spin effects in the hadronic production of top 
quark pairs at NLO in the strong coupling. 
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FIGURE CAPTIONS 

Fig. 1. Dimensionless scaling functions g'^qq{vi) (dotted), g^qq{vi) (full), and gqq'{7]) 

(dashed) that determine the expectation value {Oi)qq. 

Fig. 2. Same as Fig.l, but for (C2)gq- 

Fig. 3. Same as Fig.l, but for {Os)qq. 

Fig. 4. Same as Fig.l, but for {Oijqq. 

Fig. 5. Same as Fig.l, but for {0^)qq. 



8 




9 




Figure 2: 
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Figure 3: 
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Figure 4: 
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